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NEAR-FIELD  ANTENNA  MEASUREMENTS  ON  A  CYLINDRICAL  SURFACE: 
A  SOURCE  SCATTERING-MATRIX  FORMULATION' 

Arthur  D.  Yaghjian 

The  theory  for  probe-corrected  measurement  of  antennas  by  scanning  on  a 
circular  cylindrical  surface  enclosing  the  test  antenna  in  the  near-field  is 
formulated  from  a  source  scattering  matrix  description  of  the  test  and  probe 
antennas.     The  basic  transmission  formula  is  derived  without  recourse  to 
reciprocity,  and  from  a  common  center  approach  which  separates  as  an  isolated 
problem  the  probe  characterization  and  transformation.     Moreover,  it  is  shown 
how  an  experimental  technique  can,  in  principle,  determine  the  required 
transformed  probe  coefficients  without  the  use  of  addition  theorems.  Computer 
inversion  of  the  transmission  formula  is  accomplished  accurately  and  efficiently 
with  the  aid  of  the  sampling  theorem  and  FFT  algorithm. 

Key  words:  Cylindrical  scanning;  near-field  measurements;  source  scattering 
matrix. 

1.  INTRODUCTION 

The  ease  with  which  many  antennas  can  be  scanned  on  a  cylindrical  surface  by  moving  a 
probe  along  a  single  axis  for  different  azimuthal  orientations  of  the  test  antenna  about  a 
parallel  axis  makes  cylindrical  scanning  a  very  attractive  near-field  measurement  technique. 
In  addition  to  its  operational  simplicity,  cylindrical  near-field-far-field  deconvolution 
utilizes  the  FFT  in  a  straightforward  fashion  without  the  requirement  of  special  measure- 
ment probes,  and  the  far-field  emerges  in  conventional  elevation  and  azimuth  angles,  the 
latter  of  which  covers  a  full  360°.     And,  in  fact,  both  the  2-dimensional  and  3-dimensional 
formulation  of  the  probe-corrected  technique  for  measuring  antennas  on  a  cylindrical  surface 
have  been  derived  and  verified  through  experiment  by  Brown  and  Jull  [1]^  (2-D  case)  and 
Leach  and  Paris  [2]    (3-D  case).     In  both  papers  the  EM  fields  of  the  test  antenna  were 
expressed  in  terms  of  cylindrical  waves  and  the  Lorentz  reciprocity  theorem  was  used 
(assuming  reciprocal  probes)  to  derive  the  basic  "transmission  formula"  needed  to  apply  the 
cylindrical  scanning  technique.     Brown  and  Jull  expressed  the  receiving  characteristics  of 
the  probe  in  terms  of  a  plane-wave  spectrum,  while  Leach  and  Paris  used  cylindrical  waves 
to  express  the  receiving  characteristics  of  the  probe.     In  the  present  work,  the  basic 
transmission  formula  (25)  is  derived  from  the  "source  scattering  matrix"  equations  for 
the  test  and  probe  antennas  referred  to  a  common-center  cylindric?.!  coordinate  system.  This 
cylindrical  scattering  matrix  approach,  which  is  analogous  to  the  plane-wave  scattering 
matrix  approach  used  by  Kerns  [3],  proves  simple  and  straightforward,  and  does  not  require 
the  probe  or  test  antenna  to  be  reciprocal,  although  reciprocity  is  expressed  conveniently 
as  a  relationship  between  elements  of  the  scattering  matrix  (see  eqs.   (15  and  19)).  The 
parallels  between  the  cylindrical  source  scattering  matrix  formulation  used  throughout  these 
notes  and  the  plane-wave  scattering  matrix  formulation  introduced  by  Kerns  is  striking  and 
is  emphasized  to  facilitate  the  reader's  understanding  of  the  cylindrical  formulation. 

Another  distinct  advantage  of  the  source  scattering  matrix  approach  is  that  it 
clearly  and  completely  separates  the  problem  of  determining  the  receiving  characteristics 

'''This  Technical  Note  was  originally  written  as  a  tutorial  set  of  notes  for  the  NBS  short 
course  (July  7-11,  1975  -  Boulder,  Colorado). 
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Figures  in  brackets  Indicate  the  literature  references  at  the  end  of  this  paper. 


of  the  probe  from  the  derivation  of  the  transmission  formula.     The  characterization  of  the 
probe  can  be  dealt  with  as  an  isolated  topic.     It  is  irrelevant  to  the  basic  formulation 
whether  the  required  receiving  functions  for  the  probe  are,  e.g.,  measured  directly,  trans- 
formed from  its  known  far-field  or  plane-wave  spectrum  [1],  or  transformed  from  its  known 
cylindrical  mode  expansion  [2].     Of  course,  the  accurate  characterization  of  the  probe  is  a 

requirement  of  any  reliable  probe-corrected  theory JaadiS-^iseitssed  as  a  separate  topic  in 

•  usn'iol  0oi8ex!nansi3  ol  ■; 
section  4  below.  .     ,    ^  ... 

::'^fr''''THE  RADIATING  FIELDS  OF  AN  ANTENNA  IN  CYLINDRICAL  WAVES 

This  work  addresses  the  problem  of  measuring  the  radiating  fields  of  a  test  antenna  by 
scanning  with  a  probe  antenna  on  an  imaginary  cylinder  surrounding  the  test  antenna.  The 
problem  of  measuring  the  receiving  characteristics  of  the  test  antenna  is  very  closely 
related,  but  for  the  sake  of  brevity  it  will  be  left  as  an  exercise  for  the  reader.  Of 
course,  most  antennas  are  reciprocal  so  once  the  radiating  or  transmitting  characteristics 
are  determined  the  receiving  properties  follow  immediately  from  the  reciprocity  relationship 

,M^^if:>a:j  insms-xu^Bsm  bLsn-iBBn  3VX.:!0B-s:S:iB  yisv  s  gnlansoa  Inol-ibxiUxo  a±xs  isil&isq 

.As.  a  first  step  in -formulating,  t^he  cyJ-indrical  scanning  technique,  the  EM  fl-eild-^; -out-slide 
an  arbitrary  test  antenna  radiating  into  free  space  will  be  expanded  in  a  complete  s&tr  &fcjj 
cylindrical  waves  (eigenf unctions  or  modes)  .     The  test  antenna  _and.  its  j^pordinat^^ 
fixed  to  the  test  antenna  are  ^sh.ovm  _s.chematically  in  f  igure-pL.r  i  In  .cylojadrical:  ,G;oqr.dina„^^^^ 
the  point  P  is  designated  by_^(p_,  J)j._z)_f__  In  spherical  coordinates  by  (r,  (fi,  6).-  a.oxjslunrxo'i 
It  has  been  shown  [4]  that  in  cylindrical  coordinates  the  (E,  H)  fields  in  free  space; 
(in  this  case  outside  the  radius  (,a)  .of  the  ^smaj.^lest  cylinder  circumscribing  the:^  teiSft  antenna) 
can  be  divided  into  a  TE  and  TM  part  with  respect,  to,  the  z-direction.     Sp.e,c-if iQa^l>y£jaa3-.ra,...3 

..  .,r_,,.,,,  -j^.^po..,,,  °         .'^.-.-Tn  frit  ic  .;,:;:vx33Sa  sil^  sestJv 

where        and.  ij;2.,  which  give  rise  to  the  TE  and  TM  parts  of  the  field  respfe'Stively;^^"^"'"'^"'^^'^ 

..satisfy  the  scalar  wave  equations.:.'  -xsjnsD-nosimoo  &  oi  ha-xiaisi  z&nue2n&  adosq  bna  3as3  axii 
oni.sl....  .v.w-s.BXq  sn.       ^^ol^^  .4  . ^-^-^q.^  xi..B.  a--«Boe  I.:-.--^.^ 

All  quantities  with  hats  ,(^)  over  thgm  denote r unit,. vectors  ,  e^^'  (u)?fQ^-i  feimepdependeaceQ  sxi;! 
and  the  rationalized^ mks  system  .are  used  throughput,  and  the  free-space  wave  number  k  is  be 
defined  as  usual  by  k  =  (jj/iTT"  =  ^  =^|j  ^,3.035 jsoa  soiuos  iBol-ibnilxo  sri3  nssw^sd  alsXIs-xeq 

fj^-in  The  solution  to  the  scalar  wave  equation  (2)  in  eyliridf idai'-dodrditnateis^^^'^^  the 
antenna  1  is  written  in  terms  of  the  well-known  scalar  cylindrical  waves 

.~TT:o-5y.;i  a  as  nsrUxi^r  YXIft^JMi"o  esv  ..  .  . 


j-xorls  eivj  ~oi  8s3on 

.(obBttoIoD     isfc-  iy2 
K  =  (k^  -  Y  ) 


LuL.)  98:ri!o.-i 


s  (  \ 

The  B  (y)  are  arbitrary  constant  coefficients  and  H^I'^(kp)  are  the  Hankel  functions  of  the 
n  n 

first  kind.     The  Hankel  functions  of  the  first  kind  are  the  only  cylindrical  functions  which 

possess  the  proper  behavior  as  p-x"  for  k  chosen  positive  real  or  imaginary.     (Later,  when 

the  probe  enters  the  picture  giving  rise  to  sources  outside  p,  Bessel  functions  J  (kp)  are 

n 

also  Included  in  the  complete  set  of  cylindrical  modes.)     When  the  solutions  (3)  for  il;  are 

s 

substituted  into  eqs.   (1)  and  the  various  curls  are  taken  in  cylindrical  coordinates,  the 
(E,  H)  fields  can  be  written 


E(p,<^,z)  =     I      j  [B^(y)  +  b2(y)  N^l^dY  (4a) 


M^l^  +  b2(y)  N^l^ 
nY         n  nY 

(p,<t>,z)  =  ^    I         /   [B^(y)  N^l^  +  b2(y)  M^l^ldY  . 


o  n=— °°  — "» 


(Ab) 


(Z    =  /y  /e    "  377  ohms) 
o          o  o 


where  the  M  and  N  functions  are  given  by 


M^-^' (p,4),z)  =  V  X  (H^^'(Kp)e     ^e  '  6  )  =  [ —  H^^'(kp)§    -  k  H^^'     e^]e     ^e  '  (5a) 
nY  n  .j^yg  g,  z  p      n  p  n  ()> 

N^^^p,*,z)  =  iv  X  ^(l)  =  i[iY<  H^l)  V,-  "^.^^h  ^  ^^H^l^e  ]e^"*e^^^  .  (5b) 
nY  k"'--  -'^'^•ti'Y         k  n   i:_pji-f._P-i:ii-ja:;.-../.(ti--         n  z 

g 

If  the  B^(y)  coefficients  were  known  for  the  test  antenna,  eqs.   (4)  and  (5)  would  give  the 

electromagnetic  field  of  the  test  antenna  everywhere  outside  the  radius  of  the  smallest 

cylinder  circumscribing  the  test  antenna.     Essentially  the  main  purpose  of  this  work  is  to 

s 

formulate  a  procedure  for  finding  the  B^(y)  for  an  unknown  test  antenna,  by  scanning  on  a 
cylindrical  surface  with  an  arbitrary  but  known  probe  antenna.     The  only  quantities  in 
eqs.   (4)  which  change  from  antenna  to  antenna  are  the  coefficients  B^(y). 

"'^ 'If' one  Is  fortunate  enough  to  have  a  probe  at  his  disposal  which  approximates  an  ideal 
dipole,  i.e. ,  a  probe  which  measures  the  transverse  components  of  the  electric  field,  then 
the  B^(y)  can  be  found  directly  and  almost  trivially  from  eq.   (4a)  alone.     To  show  this,  the 
following  orthogonality  relationships  for  the  M  and  N  functions  are  required 

CO  ^TT 

f  f  (M^l^  X  M^}\)-6  d(fi  dz  =  0  (6a) 
( a  c ,)           J   '       nY         n  Y  P 

o  ' :  -  - 

:  dulBV  s    /     f  (N^l^  X  N^}\)'e  d(()  dz  =  0  (6b) 
nY         n  Y  P 

a„  -<»  o 

1    /(N^l^  X  M^V  ,)  -^  d.))  dz  =  "^^^(^P)"^^"*' (^P)  ^       ,<5(y+y')     .  (6c) 

-'''       nv         nv        p  k        n  n  -nn 

-00  o 

These  relationships  derive  from  the  definitions  of  M  and  N  in  eqs.  (5). 

Now  cross  N  into  eq.   (4a)  and  use  the  orthogonalities  (6)  to  yield  the  B^(y), 
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99S)  4Tr2K3H'^l'' (kp  )h'^1''(kp  )  -«>  o  -n,-Y 


Bm-xe.  nx  ^  ^    j  ,  E(p^, * , z)) -.^d*  d. 
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or  written  out  in  component  form, 


Att^k^H^^^  (kp  )  -° 


n 

-00  o 


b1(y)  =   =fT^   /  /  E  (p   ,(|,,z)  +  k2e  rp^,^,z)     e-^"V^^^d(t,dz     .  (7a) 


p  z  o 
o 


p^  is  the  radius  of  the  imaginary  scanning  cylinder  on  which  the  transverse  electric  field 
is  measured.     Similarly  cross  M  into  eq.    (4a)  to  yield  the  B^(y), 

Sn(Y)  =   ^1   /  Te  (p  .*,z)e-^"V^^^  d<^dz     .  (7b) 

n  o 

Equations  (7)  give  the  required  modal  coefficients  in  terras  of  the  measured  transverse 
electric  field  on  a  cylinder.    Moreover,  the  sampling  theorem  and  fast  Fourier  transform 
(see  section  5)  can  be  utilized  for  the  efficient  computation  of  the  double  integrals  in 
eqs.   (7).     Thus  for  an  ideal  dipole  probe  the  deconvolution  problem  reduces  to  a  trivial 
application  of  the  orthogonality  relationships  for  the  M  and  N    functions.     This  is  true  for 
planar  and  spherical  scanning  as  well — and  as  a  matter  of  fact,  for  any  of  the  six  separable 
coordinate  systems  for  which  M  and  N  vector  wave  solutions  exist  [4]. 

2.1.     Expressions  for  the  Far-Field,  Power  Gain 
Function,  and  Complex  Polarization  Ratio 

Equations  (4)  and  (5)  give  the  electromagnetic  fields  everywhere  outside  the  test 
antenna  radiating  into  free  space.     Often,  however,  only  the  far-fields  of  the  test  antenna 
are  of  interest.     Thus,  it  is  helpful  to  derive  explicit  expressions  for  the  far-field  from 
eqs.    (4)  and  (5).     This  has  been  done  in  appendix  A  by  a  straightforward  manipulation  of  the 
Fourier  transforms  without  using  asymptotic  methods  such  as  the  method  of  stationary  phase 
or  steepest  descent.     The  final  results  written  in  the  spherical  coordinates  of  figure  1  are: 

oo 

—  V  •  A 

E(r,(f,e)  =  -  2  k  sin  6   ^       (-i)'^[B^(k  cos  e)e,  -  iB2(k  cos  6)e^]e^^'''  (8a) 

r  n=-'=°  n  on  o 

e  ^E  ikr 
I(r,^,e)  =    I  ,  2  k  sin^e  e  ^     (-i)^[Bl(k  cos  9)6^  4-  iB^  (k  cos  e)e  ]e""*  .  (8b) 

o  o  n=-°° 

The  far-field  equations  (8)  do  not  involve  an  integration,  only  a  summation,  and  the  values  of 

S  I     I  s 

B^(y)  for   I Y I >k  do  not  affect  the  far-fields.     Note  also  that  at  least  some  of  the  B^(y)  must 

behave  as  1/sin  6  as  if  the  far-field  patterns  are  finite  there. 

The  power  gain  function  is  defined  as 

G(OA)  =  '"f  fl'^"^  ,  (9a) 
o  input 

where  P.  is  the  net  power  input  to  the  antenna.     The  input  power  can  be  written  in  terms 

of  the  ingoing  (a^)  and  outgoing  (b^)  modes  of  the  waveguide  which  feeds  the  antenna  (see 


figure  1  and  Kerns  [3]), 


P 

inp 

where        is  the  characteristic  admittance  of  the  waveguide  and        is  the  reflection 

coefficient  at  the  reference  plane  S  (b    =  r  a  ) .     Substitution  of  P  and  the  far- 

o    o        o  o  input 
electric-field  from  eq.   (8a)  into  eq.    (9a)  yields  the  final  expression  for  the  gain 

function. 


=  iVla'|^(l-|r   P)^l  ^     (-i)X(^  e}e^'"^\^  +  I   I  (-i)V(k  COS  e)e^"*|2} 


(9b) 


The  complex  polarization  ratio  defined  here  as  E^/E^  may  be  written  directly  from  eq.  (8a), 

-i    I  (-i)V(k  cos  e)e^^* 
6  n=-oo 

Y'  =  ~   •  (10) 

*        J  (-i)^Bl(k  cos"e)e^^* 
n 

2.2    Total  Power  Radiated 

Also,  from  the  far-field  equations  (8)  it  is  a  fairly  simple  matter  to  find  the  total 
power  P^  radiated  by  the  test  antenna, 

  1        2Tr  _ 

P^  =  2  ^  (ExH*)-§^da  =  2^  /  /   |E|2r2sin  6  d^dd  .  (11a) 
r-«o  o  o  o 

After  substituting  the  electric  field  from  eq,    (8a)  into  eq.   (11a)  and  making  use  of  the 
orthogonality  of  the  e"'"'^"^  functions,  P^  can  be  written 

2        CO  -j^ 

P    =  -^p-    I      /  [lBl(k  cos  e)  [2  +  |B2(k  cos  0)  |2]sin3e  d6.  (lib) 
o    n=-«'  o 

A  change  of  the  variable  of  integration  from  k  cos  9  to  y  transforms  eq.   (lib)  into  the  final 
expression  for  the  total  radiated  power, 

P,.=^I         /  K2[|B^(Y)i2  +  |B2(Y)|2]dY     .  (12) 
o  n=-<^  -k 

Observe  that  no  power  is  carried  in  the  waves  for  [yl^k. 

3.     DERIVATION  OF  THE  TRANSMISSION  FORMULA 


In  section  2  we  derived  the  E  and  H  fields  everjwhere  outside  the  test  antenna  in  terms 
of  its  cylindrical  wave  coefficients  B^(y) •     Explicit  expressions  were  bitten  for  the  far- 
fields,  power  gain  function,  complex  polarization  ratio,  and  total  power  radiated. 
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In  this  section  we  will  derive  the  basic  transmission  formula  needed  to  compute  the 
coefficients  B^(y)  for  an  unknown  test  antenna  by  scanning  on  a  cylinder  with  a  probe  of 
known  receiving  characteristics.     The  approach  that  is  taken  is  to  first  write  "source 
scattering  matrix  equations"  for  the  test  antenna  with  respect  to  the  coordinate  system 
(previously  defined  in  figure  1)  fixed  in  the  test  antenna.     Secondly,   the  source  scattering 
matrix  equations  for  the  probe  are  written  for  a  cylindrical  coordinate  system  (p',!f>',z') 
fixed  to  the  probe  but  with  its  z'-axis  lying  along  the  z-axis  of  the  test  antenna  coordinate 
system.     Thirdly,  the  relationship  is  derived  between  the  modal  coefficients  defined  in  the 
coordinate  system  of  the  probe  and  the  modal  coefficients  defined  in  the  coordinate  system 
of  the  test  antenna.     This  relationship  is  entered  into  the  source  scattering  matrix  equa- 
tions and  the  resulting  equations  are  combined,  neglecting  multiple  reflections,  to  give  the 
desired  transmission  formula.     Finally,   the  transmission  formula  is  inverted  to  yield  the 
B^(T). 

3.1.    The  Source  Scattering  Matrix  Equations  for  the  Test  Antenna 

Consider  the  cylindrical  coordinate  system  fixed  in  the  test  antenna  as  shown  in 
figure  1.     While  measurements  are  being  taken  with  the  probe  antenna,  it  will  act  as  a 
source  of  EM  fields  outside  the  test  antenna.     Thus  the  cylindrical  mode  expansion  in  the 
coordinate  system  of  the  test  antenna  for  radii  between  the  test  antenna  and  probe  must 
include  two  linearly  independent  radial  solutions — not  just  the  H^^^(kp)  solutions  as  was 
the  case  for  the  test  antenna  in  free  space.     As  we  shall  see,  the  present  derivation  is 
greatly  facilitated  by  choosing  J^(kp)  as  the  second  linearly  independent  solution. 

The  E  and  H  fields  between  the  test  and  probe  antennas  may  be  written  as  in  eqs.  (4), 
but  with  the  J^(kp)  solutions  included  as  well.  Specifically, 

oo  00 

E(p,<^,z)  =     I      /{[b^(Y)M^l^  +  b2(Y)N^l^]  +  [al(Y)M      +  a2(Y)N     ]}dY  (13a) 
n.=— oo  —00  IT  11 

H(p,(}),z)  =  VxE     ,  (13b) 

o 

where  M^-^\  N^^^  have  the  first  Hankel  functions  H*"^^(kp)  for  their  radial  dependence,  and 
_       _  nY        nY  n 

M    ,  N      have  the  Bessel  functions  J  (kp)  for  their  radial  dependence.     The  modes  with 
nY      nY  n 

coefficients  b  (y)  represent  the  fields  from  sources  (induced  as  well  as  applied)  inside 

s 

or  on  the  surface  of  the  test  antenna.     The  modes  with  coefficients  ^^(j)  represent  the 
fields  from  sources  (induced  as  well  as  applied)  inside  the  probe  antenna,  i.e.,  outside 
the  test  antenna. 

Also,  the  waveguide  modes  with  coefficients  a    and  b    represent  fields  from  sources 

o  o 

inside  and  outside  the  generator  (assumed  shielded)  of  the  test  antenna,  or  equivalently, 
outside  and  inside  the  test  antenna  proper. 
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In  brief  then,  the  modal  coefficients  (a^,  a^Cy))  and  (b^,  b^Cy))  represent  fields  from 
sources  outside  and  inside  the  test  antenna  respectively.     The  inside  coefficients 

S  Q 

(b^,  t>^(y))  are  related  to  the  outside  coefficients  (a^,  a^(y))  by  a  linear  matrix  trans- 
formation which  will  be  called  the  source  scattering  matrix  for  the  test  antenna,  to 

distinguish  it  from  the  scattering  matrix  for  the  H^^\  H^^)  representation  of  the  fields 
/  \  n  n 

instead  of  the  H^^,        representation  used  here.     Specifically,  these  source  scattering 
matrix  equations  are  written 

2  CO  CO 

^0  =  ^0^0+^      ^        /    ^n^^)^^^)  '^^  (1^^) 

S=l   n=-a>  -00 

2  CO  00 

b^(y)  =  T^(y)a^  +11        /  S^^(y,e)a^(B)dB     .  (lAb) 
q=l  m=-°°  -«= 

s  s  so 

The  functions  R^(y) ,  T^(y)  and  S^jJj(y>S)  embody  the  receiving,  transmitting,  and 
scattering  properties  of  the  test  antenna.  is  the  waveguide  reflection  coefficient 

defined  at  the  surface  S^.     For  a  given  cylindrical  coordinate  system  fixed  in  the  test 
antenna,  these  four  quantities  depend  only  on  the  character  of  the  test  antenna.     In  Kern's 
notation  [3]  S^^,  S^^,  S^^,  and  S^^  correspond  to  R,  T,  S,  and  1^,  respectively. 

The  source  scattering  matrix  equations  (14)  can  be  taken  as  a  definition  of  linearity 
for  the  test  antenna  or  they  can  be  derived  from  Maxwell's  equations  through  a  "source 
adjoint  reciprocity  lemma."    The  derivation  of  this  source  adjoint  lemma  and  the  scattering 
matrix  equations  (14)  may  be  found  in  appendix  B.     It  is  also  proven  in  appendix  B  that  the 
receiving  and  transmitting  functions  of  a  reciprocal  test  antenna  (or  between  a  nonreciprocal 
antenna  and  its  adjoint)  satisfy  the  relationship 

R%(-y)  =  (-1)"  ^    T^y)   ,  (15a) 
-n  n  Z  K  n 

o  o 

and  the  scattering  function  satisfies  the  relationship 

S"!^       (-B,-y)  =  S^l  (y,B).  (15b) 
-m , -n  nm 

When  the  probe  antenna  is  removed  so  that  the  test  antenna  is  radiating  into  free  space, 
eq.   (14b)  becomes 

b^(y)  =  B^(y)  =  T^(Y)a^.  (16) 

f  s) 

Thus  if  we  can  determine  T      (y)  for  the  test  antenna  we  can  compute  the  free-space  EM  fields 

n 

of  the  test  antenna  by  means  of  eqs.    (4).     And,  in  particular,  the  far-fields,  power  gain 
function,  complex  polarization  ratio,  and  total  power  radiated  could  be  computed  from 
eqs.   (8),   (9b),   (10),  and  (12),  respectively.     Thus  the  problem  of  measuring  the  radiating 
properties  of  the  test  antenna  reduces  to  that  of  determining  T^(y)-     The  next  step  toward 
this  end  is  to  write  the  source  scattering  matrix  equations  for  the  probe  antenna. 
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3.2.     The  Source  Scattering  Matrix  Equations  for  the  Probe  Antenna 


Consider  a  cylindrical  coordinate  system  fixed  in  the  probe  antenna  but  with  its 
z ' -axis  lying  outside  the  probe  antenna  along  the  z-axis  of  the  cylindrical  coordinate 

system  fixed  in  the  test  antenna  (as  shown  in  figure  2).     a'  and  b'  are  the  coefficients  of 

—d.   o  o 

the  ingoing  and  outgoing  waveguide  mode  of  the  probe  antenna,  and,  in  particular,   the  primes 
on  the  cylindrical  coordinate  system  emphasize  the  fact  that  this  coordinate  system  is  fixed 
in  the  probe.     Of  course,  the  probe  coordinate  system  is  nothing  more  than  the  test  coordinate 
system  rotated  about  the  z-axis  through  an  angle  if ^  and  translated  along  the  z-axis  through 
a  displacement  z^.     (t))^  and  z^  vary  while  scanning.) 

The  EM  fields  on  a  cylinder  between  the  test  and  probe  antennas  can  be  written  in  terms 
of  the  prime  system  fixed  in  the  probe: 

00  00 

E'(p',(j,',z')  =     y         f{[a'l(Y)M*-l^  +  a'2(Y)N^l^]  +  [b'l(Y)M      +b'2(Y)N     ]  }dY  (17a) 
'       Xi         v^y  n          nY  n         uy        n  nY 

n=-oo  -00  '  .  ' 

H' (p' ,z')  =  V'xE'.  (17b) 

itoy 

o 

These  equations  are  identical  to  eqs.    (13)  except  for  the  primed  coordinates  (p',())',z')  and 
primed  coefficients  (a',b')  replacing  the  unprimed.     The  a'  and  b'  coefficients  have 
deliberately  replaced  the  b  and  a  coefficients  respectively,  so  that  the  a'  and  b'  coeffi- 
cients would  now  be  associated  with  sources  outside  and  inside  the  probe  antenna 
respectively.     In  other  words,  the  a'  now  refers  to  the  H^^^(kp')  modes,  and  the  b'   to  the 

J  (kp')  modes, 
n 

Just  as  was  done  for  the  test  antenna,  the  inside  probe  coefficients  (b' ,b'  (y))  can  be 
related  to  the  outside  probe  coefficients  (a^,  a^  (y))  by  a  source  scattering  matrix  for  the 
probe : 

K  =  ^Yo  +    I      I  a;^Y)dY  (18a) 

s-1  n=-°°  -<» 

2  CO  CO 

b'^Y)  =  T'^a'  +    I      I        /s'^^(Y,e)a'^  dB     .  (18b) 
n  no  ^        ''nm  m 

q=l  n=-°o  -oo 

s  s  so 

The  functions  (y) ,  T^  (y) ,  and  ^^j^  vY»^)  represent  receiving,  transmitting,  and  scattering 
properties  of  the  probe  antenna,  is  merely  the  reflection  coefficient  for  the  probe 

waveguide  referred  to  the  surface  S\     For  a  given  cylindrical  coordinate  system  fixed  in 
the  probe  antenna,  these  four  quantities  depend  only  on  the  character  of  the  probe  antenna. 

Like  eqs.   (14),  the  probe  eqs.    (18)  are  derived  from  Maxwell's  equations  in  appendix  B. 
Appendix  B  also  proves  that  the  source  scattering  matrix  for  a  reciprocal  probe  (or  for  a 
non-reciprocal  probe  and  its  adjoint)  obeys  the  following  relationships  analogous  to 
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eqs.   (15)  which  apply  to  the  test  antenna, 


o  o 

The  probe  scattering  equations  (18)  can  be  linked  directly  to  the  corresponding  test  antenna 

scattering  equations  (14)  by  relating  the  coefficients  (a'^,  b'^)  of  the  cylindrical  modes 

s  s 

in  the  coordinate  system  fixed  in  the  probe  to  the  coefficients  (a^,  b^)  of  the  cylindrical 
modes  in  the  coordinate  system  fixed  in  the  test  antenna.     This  relationship  between 
(a^^,  b^^)  and  (a^,  b^)  is  surprisingly  simple  and  can  be  found  in  a  surprisingly  simple 
fashion. 


(19a) 
(19b) 


3.3.     Relationship  Between  the  Probe  and  Test  Antenna 
Coefficients  in  the  Common-Center  Systems 

We  know  that,  regardless  of  whether  the  probe  or  test  antenna  coordinate  system  is  used, 
the  EM  fields  at  every  point  in  space  must  be  the  same,  i.e.. 


E'(p',  (}.•,  z')  =  E(p,  (j>,  z)  (20a) 

H'(p',  (|)',  z')  =  H(p,  (]),  z)   ,  (20b) 

when  (p',  (})',  z')  and  (p,  ((),  z)  refer  to  the  same  point  in  space.     But  as  the  probe  scans  the 
test  antenna,  the  probe  coordinate  system  is  merely  rotated  through  an  angle  tj)^  about  the 
z-axis  and  translated  a  distance  z^  along  the  z-axis,  i.e., 

p'=p,  (()'=<})  -  i})^,  and  z'  =  z  -  z^ 

and  eqs.   (20)  become 

E'(p,  4.  -  (}>^,  z  -  z^)  =  E(p,  4),  z)  (21a) 

H'(p,  (f)  -  z  -  z^)  =  H(p,  4,  z)   .  (21b) 

Now  when  (E,  H)  and  (E' ,  H' )  are  written  explicitly  in  terms  of  the  linearly  independent 
cylindrical  waves  in  eqs.   (13)  and  (17),  the  only  way  eqs.    (21)  can  be  satisfied  is  if 

in(})  iyz 

a;^^(Y)  =  b^(Y)  e      °  e      °  (22a) 
in(()  iyz 

b^^(Y)  =  a^(Y)  e      °  e      °  .  (22b) 

These  are  the  desired  relationships  linking  the  probe  and  test  antenna  scattering 
equations  (18)  and  (14). 
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3.4.     The  Transmission  Formula 


Equations  (22a)  recast  the  probe  receiving  equation  (18a)  into  the  form 

2      0°        °°  int})  iyz 

^o^'^o'^)  =  ^o^o        ^      ^      /     R;"(Y)b^(Y)  e      °e      °  dy.  (23) 
8=1  n=-«'  -°° 

Moreover,  ti^(Y)  from  eq.    (14b)  can  be  substitued  in  this  eq.    (23)  to  yield 

2      00        OS  2      °o        <»  in4)  iyz 

s=l  n=-<»  -°°  q=l  m=-°°  -°° 

(24) 

S  Q 

The  contribution  to  ^^('t'Qj  2^)  from  the  term  in  eq.    (24)  containing  ^j^j^J^YjB)  represents 
multiple  reflections.     That  is,  it  is  that  part  of  the  output  of  the  probe  caused  by  fields 
which  have  left  the  test  antenna,  reflected  from  the  probe,  re-reflected  from  the  test 
antenna,  and  returned  to  excite  the  probe.     This  process  repeats  itself  ad  infinitum. 
Assuming  these  multiple  reflections  are  negligible ,  and  the  input  a^  to  the  probe  is  kept 
zero,  eq.    (24)  reduces  to 

CO    2  intf)  iyz 

^o^^o'^o^  =  %     ^  ^  R;"(Y)T^(Y)]e      °  e      °  dy     .  (25) 

n=-oo  -00  s=l 

(If  a^  is  not  kept  zero  but  the  probe  is  terminated  in  a  load  with  reflection  coefficient 

r'     the  right  side  of  eq.    (25)  is  merely  divided  by  (l-rT').) 

(1)       (2)  °  (1) 

Here,  it  can  be  noted  that  if  (H       ,  H      )  functions  had  been  used  instead  of  (H      ,  J  ) 

n          n  n  n 

for  the  expansion  of  the  fields  in  cylindrical  waves,   it  would  be  impossible  without  going 
through  an  extremely  tedious  procedure  to  show  the  effect  of  neglecting  multiple  reflections 
since  H^"*^^  waves  would  be  excited  by  the  probe  as  well  as  test  antenna. 

This  "transmission  formula"  (25)  relating  the  output  of  the  probe  and  the  "coupling 

2 

product"  (J    R'    (y)  T  (y)  where  R'    (y)  is  the  receiving  function  of  the  probe  and  T  (y) 
s=l 

is  the  transmitting  function  of  the  test  antenna)  is  analogous  to  the  "transmission  integral" 
of  the  planar  formulation  [3],  and  forms  the  central  relationship  on  which  the  cylindrical 
scanning  technique  is  based.     In  essence,  it  is  the  same  formula  as  eq .    (23)  in  Leach  and 
Paris  [2],  and  the  same  formula  as  the  2-dimensional  equation  (18)  in  Brown  and  Jull   [1]  if 
the  polarization  index  s,        -  dependence,  and  y  integration  are  deleted  from  eq.  (25). 

The  transmission  formula  (25)  has  been  derived  in  a  simple,  straightforward  manner  from 
the  source  scattering  matrix  equations  for  the  test  and  probe  antennas  under  the  assumption 
of  negligible  multiple  reflections.     No  other  restrictive  assumptions  were  made      the  probe 
and  test  antenna  may  even  be  non-recrlprocal .     (The  derivations  of  Leach  and  Paris   [2]  and 
Brown  and  Jull  [1]  inherently  assume  a  reciprocal  probe  although  both  their  derivations  can 
be  modified  slightly,  using  adjoint  reciprocity,   to  apply  to  non-reciprocal  probes  as  well.) 


10 


3.5.     Inversion  of  the  Transmission  Formula 


Let's  take  a  closer  look  at  the  transmission  formula  (25).     In  principle,  as  the  probe 

scans  on  the  surface  of  an  imaginary  cylinder  surrounding  the  test  antenna,  the  output 

b'((i>  ,z  )  of  the  probe  is  recorded  for 
o    o  o 

0  <  (fi     <  Zir 
—  o 

-»  <  z     <  <» 
o 

That  is,  the  amplitude  and  phase  of  b' (A  ,  z  )  is  measured  for  all  values  of  6    and  z  .  (In 

o     o       O  0  0 

practice,   the  z^  scan  is  limited  to  some  finite  scan  length  and  b^((t>^,  z^)  is  assumed 

negligible  outside  this  region.     Also,  data  need  be  sampled  and  recorded  only  at  a  finite 

number  of  measurement  points.     These  topics  of  data  collection  and  computations  will  be 

discussed  in  section  5  below.)     Having  measured  b^((t)^,  z^)  >  the  Fourier  series  and  integral 

of  eq.   (25)  can  immediately  be  inverted  to  yield  the  solution  for  the  coupling  product  in  terms 

of  the  measured  data  b'((j)  ,  z  ): 

o    o  o 

2  «=  2Tr  -in())  -iyz 

1  R;^y)  <(Y)  =  4^    /  /  b;(*,,z^)e        °  e        °  d*^dz^     .  (26a) 
s=l  O  o 

By  computing  the  double  integral  on  the  right  side  of  eq.    (26a),  the  value  of  the  coupling 

product, 

r;^(y)t1(y)  +  r;^2(y)x2(y)  , 

can  be  determined  for  as  many  n  and  y  as  desired.     However,  even  if  the  receiving  charac- 
teristics (R' ^ ,  R' 2)  of  the  probe  are  known,  the  transmitting  characteristics  (T^,  T^)  of 
n        n  ^  n  n 

the  test  antenna  are  not  determined  uniquely  (2  unknowns  for  each  equation) .  Fortunately, 
all  that  is  needed  is  another  scan  with  a  different  probe,  or  the  same  probe  in  a  different 
orientation  (say  rotated  about  its  baseline  axis  by  90°).     In  other  words,   two  linearly 
independent  scans  are  necessary  to  account  for  the  polarization  of  the  EM  fields.  The 
second  scan  produces  a  second  equation  to  complement  eq.    (26a).  Specifically, 

2  ^  2-n  -incf)  -iyz 

y  R"^(y)T^(y)  =  \  !  b"(*  ,z  )e        °  e        °  d*    dz       ,  (26b) 

'-^  n     '     n  '        AiT'^a      ^   '     o    o    o  o  o 

s=l  o  -«=  o 

where  R"^(y)  and  b"(d),  z  )  are  the  receiving  characteristics  and  output  of  the  second  or 
n  o  o 

reoriented  probe  respectively. 

Assuming  for  the  moment  that  the  receiving  functions  (R' ,  R")  of  the  probe(s)  are  known, 
eqs.   (26)  can  be  solved  immediately  for  the  transmitting  functions  T^(y)  of  the  test 
antenna : 

t1(y)  =  [R;;2(Y)r  (y)  -  r;2(y)i;;(y)]/a^(y)  (27a) 

t2(y)  =   [R'2(y)I"(y)  -  R"1(y)I'(y)]/A^(y),  (27b) 
n  n         n  n         n  n 
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where  I' (y)  and  I"(y)  are  the  integrals  on  the  right  side  of  eqs.    (26a)  and  (26b) 


n  n 


respectively;  and  A  (y)  is  the  determinant 


=  R"2(y)R'1(y)  -  R'2( 

n         n  n  n 


(28) 


Provided  the  receiving  characteristics  of  the  probe(s)  are  known,  eqs .    (26) ,   (27) ,  and  (28) 


Naturally  the  question  arises  as  to  how  the  receiving  characteristics  of  the  probe  can 

s  s 

be  determined.     Specifically,  how  do  we  find  the  R'   (y)  and  R"  (y)  which  are  required  in 

s  s 

eqs.    (27)  and  (28).     (In  the  subsequent  discussion,  I  will  mention  only  since  R^  can 

be  found  by  the  same  procedure.) 


transformations  but  requires  two  cylindrical  scans  and  two  test  antennas  (or  one  test 
antenna  in  two  orientations)  with  known  far-fields.     The  second  method,  which  is  closely 
related  to  the  approach  used  by  Brown  and  Jull  [1]  and  Leach  and  Paris   [2]  for  cylindrical 
scanning,  and  by  Jensen  [5]  and  Wacker  [6]  for  spherical  scanning,  involves  a  coordinate- 
like transformation  between  the  known  receiving  functions  of  a  probe  in  a  coordinate 
system  centered  on  the  probe  and  the  required  receiving  characteristics  R^  (y)  which  have 
been  defined  with  respect  to  a  coordinate  system  centered  outside  the  probe. 

The  chief  advantage  of  the  first  method  is  that  it  requires  no  coordinate  transforma- 
tions, which  can  be  rather  cumbersome  expressions.     Its  chief  disadvantage,  besides  the 
fact  that  it  requires  the  a  priori  knowledge  of  the  far-fields  of  a  test  antenna,  may  be 
that  the  given  antenna  must  have  non-negligible  transmitting  functions  T^(y)  foi^  all  values 
of  n  and  y  required  of  the  probe.     This  will  be  explained  specifically  below.     The  second 
method  consisting  of  transforming  the  receiving  functions  of  a  known  probe  does  not  exhibit 
this  limitation.     It  should  be  mentioned  that  the  two-identical-antenna  or  generalized- 
tliree-antenna  techniques  sometimes  applicable  in  planar  near-field  scanning  [3]  do  not 
appear  feasible  for  cylindrical  scanning  because  of  the  more  complicated  coordinate  trans- 
formations (such  as  eqs.    (37)  and  (38)  below)  which  are  involved  in  the  cylindrical  case. 

In  any  case,  it  is  emphasized  that  the  common-center,  source  scattering  matrix  approach 
applied  in  the  present  paper  to  the  problem  of  cylindrical  scanning  results  in  a  trans- 
mission formula  (25)  which  embodies  the  essence  of  non-planar  scanning  in  a  simple  form 
(analogous  to  the  planar  transmission  integral)  whose  derivation  was  accomplished  without 
the  introduction  of  cumbersome  coordinate  transformations.     Questions  related  to  the 
characterization  of  the  measuring  probe  and  the  associated  transformations  can  be  concen- 
trated on  and  discussed  as  a  separate  topic. 


4. 


THE  RECEIVING  CHARACTERISTICS  OF  THE  PROBE  ANTENNA 


Two  methods  will  be  explained  for  determining  R*    .     The  first  involves  no  coordinate 
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4.1.     Measuring  the  Probe  Receiving  Functions 


Using  a  Known  Test  Antenna 


This  first  method  for  determining        (y)  is  quite  easy  to  explain.     The  unknown  probe 

whose  receiving  characteristics  are  wanted  scans  two  test  antennas  (or  one  test  antenna  in 

two  orientations)  whose  far-fields  (amplitude,  phase,  and  polarization),  say,  are  known. 

Since  the  far-fields  are  given,  eqs.    (4)   (with  the  asymptotic  form  of  the  Hankel  function 

for  large  p  inserted)  or  eqs,    (8)  for   |y|  £  k  can  be  inverted  using  the  orthogonality 

relationships  (6)  between  M  and  N  to  yield  the  transmitting  functions  T^(y)  of  the  known 

s  ^ 
test  antenna(s).     With  the  T    of  the  test  antenna(s)  known,  eqs.    (27)  with  the  T's  and  R's 

interchanged  can  be  used  to  determine  the  receiving  functions  R^    of  the  probe.     To  use 

eqs.    (27),  however,  the  determinant  ^^(y)  of  eq-    (28),  now  containing  the  T's,  should  not  be 

too  small.     This  implies,  as  mentioned  above,   that  the  T's  of  the  known  antenna(s)  must  not 

be  negligible  for  a  given  n  and  y. 


4.2.     Transformation  of  the  Receiving  Functions  of  a  Known  Probe 

f  s 

This  second  method  for  determining  R^  (y)  assumes  a  probe  of  known  receiving  character- 
istics (or  transmitting  characteristics,  if  reciprocal)  but  known  in  a  coordinate  system 
centered  at  the  probe,  whereas  R^  (y)  was  defined  with  respect  to  a  cylindrical  coordinate 
system  fixed  to  the  probe  but  centered  outside  the  probe.     Here  it  will  be  assumed  that  the 
known  receiving  characteristics  are  given  in  terms  of  a  cylindrical  wave  expansion,  in  order 

to  conform  with  the  work  of  Leach  and  Paris   [2].     The  transformation  to  R'    (y)  from  known 

n 

receiving  characteristics  given  in  terms  of  a  plane  wave  spectrum  has  been  derived  by  the 
author  but  the  derivation  will  not  be  included  in  these  notes.     This  latter  transformation 
is  similar  to  that  derived  by  Brown  and  Jull  [1]  for  the  2-dimensional  cylindrical  problem, 
and  recently  by  Borgiotti  [11]  for  the  3-D  problem.     As  Borgiotti  points  out,  one  advantage 
of  the  plane-wave  transformation  is  that  the  integrals  involved  lend  themselves  to  asymptotic 
evaluation. 

g 

Consider  a  probe  whose  receiving  functions  ^^^^(y)  ^^^^  known  with  respect  to  a  cylindri- 

g 

cal  coordinate  system  fixed  in  the  probe  and  centered  on  itself  (as  opposed  to  R^  (y)  which 
is  defined  with  respect  to  a  cylindrical  coordinate  system  fixed  in  the  probe  but  centered 
outside  the  probe).     Let  the  z-coordinates  (but,  of  course,  not  the  z-axes  which  are  separated 
by  a  distance  d)  be  the  same  for  each  system.     Also  choose  the  baseline  from  which  the 
angles  in  each  system  are  measured  to  coincide.     The  coordinate  system  centered  on  the  probe 
will  be  called  the        system  and  will  have  cylindrical  coordinates  labelled  (p^^,  (j)^,  z^) . 
The  coordinate  system  centered  outside  the  probe  at  the  test  antenna  but  still  fixed  in  the 
probe  will  be  called  the  C  system, and  as  before,  its  coordinates  will  be  labelled 
(p' ,())' ,z')  .     Refer  to  figure  3  for  a  schematic  of  the  two  systems. 
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The  relationship  between  the  output  of  the  probe  (for  a'  =0)  and  the  receiving 
functions  ^^-^n^^^  ^1  ^^^^^"^  definition  (compare  with  eqs .  (14)), 

2      °°  ^ 

=     J       J       /  R;^(Y)a;^(Y)dY  ,  (29) 
s=l  n=-oo  -oo 

where  a'     are  the  modal  coefficients  of  the  J  (kp^)  modes  which  are  excited  by  sources 
in  n  i 

existing  outside  the  probes,  i.e.,  by  sources  on  the  test  antenna.     The  values  of  these 

g 

^In           ^^i^'^t^  independent  of  sources  within  p^,  i.e.,  sources  on  the  probe. 

The  output  of  the  probe  b'   for  a'  =  0  is  given  alternatively  by  eq.    (18a)  in  terms  of 

the  required  probe  receiving  functions        (y)  in  the  C  system, 

2  00  00 

b;  =    ^      );        /  R;^Y)a;^Y)dY  .  (30) 
s=l  n=-«>  -0° 

A  comparison  of  eqs.    (30)  and  (29)  shows  that  the  required  R'    (y)  can  be  found  in  terms  of 

s  s  s 

the  given  R'    (y)  if  the  a'     can  be  written  in  terms  of  the  a'    ,  i.e.,   the  J  (kp,)  modal 
In  in  n  n  i 

coefficients  in  terms  of  the  (kp')  modal  coefficients.     Fortunately,  this  can  be 

accomplished  in  a  straightforward  manner  by  equating  the  EM  fields  written  in  the  and 
C  coordinate  systems. 

In  the        coordinate  system,   the  EM  fields  can  be  expanded  in  its  complete  set  of 
cylindrical  modes.     Likewise,  in  the  C  coordinate  system  the  same  EM  fields  can  be  expanded 
in  its  complete  set  of  cylindrical  modes.     The  fields,  regardless  of  which  expansion  is 
used,  are  the  same.     So  the  E  and  H  fields  written  in  each  coordinate  system  can  be  equated. 
Before  doing  this,  however,  a  subtle  point  should  be  explained. 

The  C  cylindrical  mode  expansion  is  valid  only  up  to  a  radius  p '   touching  the  probe 
antenna,  and  the  C.^  cylindrical  expansion  is  valid  only  up  to  a  radius  p^  touching  the  test 
antenna.     Thus,  at  first  sight,  we  might  conclude  that  the  EM  fields  cannot  be  equated  on 

any  complete  circle  of  either  the  C-j^  or  C  system.     Fortunately,   the  problem  disappears 

s  s 

under  closer  scrutiny.     Since  both  the  a'     and  a'     refer  to  modes  whose  sources  lie  at  the 

n  In 

test  antenna,  the  part  of  the  EM  fields  expressed  by  both  expansions  are  valid  for  all 

(p  ,  (}>  ,  z^)  within  any  cylinders  of  radius  p^  up  to  the  radius  which  touches  the  test 

(1)  (2) 

antenna.     Again,  it  is  noted  that  this  would  not  be  true  if  (H      ,  H      )  functions  had  been 

(1)  ^  " 

used  instead  of  (H      ,  J  )  and  the  derivation  would  have  become  unnecessarily  complicated, 
n          n  J        f  > 

if  performed  correctly. 

g 

Specifically,   the  a^    field  of  system  C  is  given  from  eq.    (17a)  as 

CO  00 

E'(p'.  z')  =    I      j     [a;\Y)M^J^p',  z')  +  a;^(Y)N^^^p',  z')]dY  •  (31a) 

n=— CO  —00 

Similarly  the  a|^  field  of  system        is  given  by 

00  CO 

E^(p^,  z^)  =     I      j     [a[^(Y)M     (p^,   c}.^,  z')  +  ain^Y)N     (P;L'  *1'  ^'^^'^^  '  ^^^^^ 

n=-°=  -o"  ' 


(z^  =  z') 
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Now  M^^^  is  defined  in  eqs.   (5)  as 
ny  -1       \  / 


m(^)(p',  z')  =  Vx  6^  Hy>(Kp')e^"*'e^^^'     .  (32) 


From  Graf's  addition  theorem  ([7],  formula  (9.1,79))  for  cylindrical  waves,  we  have  the  identity 
H^l>(Kp')e-*'  =     I      eI^I  (Kd)  JjKp^)e-*l   .  (33) 


which  allows  M^^^  in  eq.   (32)  to  be  written 
ny  n     V  / 


m(1)(p',  z')  =    J      H^^^  (Kd)M^^(p^,  <^^,   z^)     ,  (34a) 


in=-oo 


since  V'x  =  y^x  and  z'  =  z.     (That  V'x  =  v^x  follows  from  the  fact  that  the  curl  operator 
by  its  integral  definition  remains  invariant  under  coordinate  transformations.)  Similarly 

N^^^P',  z')  =     I    n^ll  (^d)V(p^,  z^)     .  (34b) 

m=-<» 

Substitution  of  eqs.   (34)  into  eq.   (31a)  yields  through  comparison  with  eq .   (31b)  (after 

interchanging  m  and  n,  and  using  orthogonality  (6))  the  relationship  between  the  J  (kp^)  and 

(1)  ^ 
H      (kp')  modal  coefficients, 
m 


,  y      H^-"^   (Kd)a'^  .  (35) 

In        ^        m-n  m  ^  ^ 

m=-«= 


Further  substitution  of  a|^  of  eq.   (35)  into  eq.   (29)  yields  (after  m  and  n  are  again 
interchanged) , 

b;  =     I      i      fll    <(Y)Hi'i(<d)]a;^y).  (36) 
s=l  n=-a>  -<»  m=-°° 

Comparing  eq.   (36)  with  eq.   (30)  while  remembering  that  the  a'     are  linearly  independent 
allows  us  to  write  the  final  relationship  between  the  desired  receiving  functions  in 
the  cylindrical  coordinate  system  centered  outside  the  probe  and  the  known  receiving 
functions         in  the  cylindrical  coordinate  system  centered  on  the  probe: 

oo 

K'^^^  =  I  •  (37) 

in=-oo 

Equation  (37)  displays  a  very  interesting  and  useful  property.     Even  though  the  probe 

s 

receiving  functions  R'     in  the  system  centered  on  itself  may  be  non-negligible  for  only  a 
Im 

few  values  of  m,  eq.   (37)  will  give  the  receiving  coefficients  R^  (y)  for  much  larger 
values  of  n.     This  means  that  the  transformation  expressed  by  eq.   (37)  can  be  applied  even 
to  very  small  probes  which  are  described  by  a  small  number  of  cylindrical  modes.  Also 
note  from  eq.   (37)  that  as  d  gets  very  small  R^  (y)  may  get  very  large. 
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Finally,  it  should  be  mentioned  that  most  probes  used  in  practice  will  be  reciprocal 
(see  eq.   (15a))  so  that  eq.    (37)  may  be  rewritten 

CO 

„,s,  .       4ttk^      V       /  TN^rpiS     ,     .„(1),  J, 

\          =  ^Tzk    I  ^l,-J-^^\-ra^^^^   '  ,  (38) 

o  o    m=— °° 

where  the  T^^(y)  times  a^  are  merely  the  modal  coefficients  of  a  probe-centered  cylindrical 

s  s 

wave  expansion  of  the  radiated  field  of  the  probe.     Of  course,  the  T'  a'  =  B'     can  be  com- 

Im  o  my 

puted  from  the  far-electric-field  by  inverting  eq.   (4a)  or  eq.    (8a)  for   |y|  using  the 

orthogonality  (6)  of  the  M  and  N  functions. 

g 

Note  that  if  the  probe  is  omni  in  the  i>  direction  T'     =  0  for  m  ^  0,  and 
s  (1)  ^  Im 

(y)  =  f(Y)H^     (i<d)  in  eq.    (38),  which  is  compatible  with  the  ideal  dipole  probe  case, 

eqs.    (7).     (f(Y)  denotes  an  arbitrary  function  of  y.) 

5.     APPLICATIONS  OF  THE  SAMPLING  THEOREM  AND  FAST  FOURIER  TRANSFORM 

Let's  assume  the  receiving  characteristics  of  the  probe (s)  have  been  determined  by 
either  the  "direct  measurement  method"  of  section  4.1.,  or  the  "transformation  method" 
using  eq.   (37)  or  eq.   (38).     Then  after  scanning  the  test  antenna  with  two  linearly  inde- 
pendent probes  (or  one  probe  in  two  orientations),  eqs.   (27)  yield  immediately  the  trans- 

s  s 
mitting  functions  T  (y)  of  the  test  antenna.     From  the  values  of  T  ,  the  far-field,  for 
o  n  '  n 

example,  of  the  test  antenna  can  be  computed  by  summing  eqs.    (8)  with 


B®(k  cos  6)  =  T^(k  cos  6)  a     .  (39) 
n  n  o  v-"^/ 

There  still  remains  an  important  question  to  be  answered.     How  closely  must  the  scan 

points  be  spaced,  at  which  the  amplitude  and  phase  of  the  probe  outputs  b^((j)^,  z^)  and 

b"((j)  ,  z  )  are  measured,  in  order  for  the  integrals  I'(y)  and  I"(y)  in  eqs.   (27)  to  be 
o    o      o  n  o 

evaluated  accurately.     This  question  is  answered  by  looking  at  the  definition  of  I^(y) 

(or  I^(y))  given  in  eqs.  (26): 

„  z  max  „  J   J.  ■ 

2  o        2tt  -in(t>  -lyz 

I'(y)  =    I  R'^(y)T^(y)  =  /  /  b'(^  ,z  )e        °e        °  dcf.  dz       .  (40) 

n  n          n  ,9  ^  .  ■'     o    o    o  00  ' 

s=l  4Tr^a    z  mm  o 

o  o 

(The  ±°=  limits  of  integration  for  z^  have  been  replaced  by  (z™"""^,  z™^^)  because  in  practice 

the  probe  output  b'(())  ,  z  )  must  be  assumed  negligible  outside  a  certain  finite  scan  length.) 

Consider  the  transmitting  functions  Tj^(y)  of  the  test  antenna.     They  represent  the 

coefficients  of  the  cylindrical  waves  emanating  from  the  test  antenna  (with  a    =  1) 

th  ^ 

transmitting  into  free  space.     The  subscript  n  refers  to  the  n      order  Hankel  function, 

H<^^^Kp)  k2  =  k2  -  Y^. 

n  ' 

It  can  be  shown  [8]  that  the  reactive  energy  in  the  cylindrical  modes  M^^\  N^^^  containing 

(1)  ny  ny 

the  Hankel  function  H^     (kp)  grows  extraordinarily  large  whenever  n^  becomes  much  greater 

16 


than  (kp)2.     This  result  implies  that  the  reactive  fields  or,  equivalently ,   the  Q  of  the 
antenna  is  extremely  large  unless  the  T^(y)  grow  extremely  small  for  n^>(,Ka)^,  where  "a"  is 
the  radius  of  the  smallest  cylinder  circumscribing  the  test  antenna .     Ordinarily,  antennas 

g 

do  not  have  a  high  Q,  and  thus  the  T^(y)  will  become  negligible  not  far  outside  the  circle 
shown  in  figure  4  and  defined  by 

n2  =  (Ka)2  =  (k2  -  y^)a^  ,  (41) 

or 

n2  +  (Ya)2  =  (ka)2  . 

There  is  a  reasonable  physical  explanation  accompanying  the  result  (41) .     Equations  (8) 

show  that  the  values  of  y  nearer  to  k  correspond  to  the  far-field  on  a  circle  (oj<<J)<2it) 

-1  Y 

nearer  to  the  z-axis  (6  =  cos      ■^)  .     Assuming  the  number  of  far-field  fluctuations  per  solid 
angle  are  roughly  the  same  in  all  directions,  the  far-field  variation  with  (J)  will  be  lesser 
for  circles  nearer  the  z-axis  (larger  y) •     Thus  the  number  of  cylindrical  modes  needed  to 
expand  the  field  will  grow  less  and  less  as  y  approaches  k,  since  the  higher  order  modes, 
which  represent  the  more  rapidly  varying  parts  of  the  field,  will  not  be  needed  to  expand 
the  slower  variations  with  <^ .     Equation  (41)  merely  expresses  quantitatively  this  intuited 
result.     Specifically,  eq.   (41)  says  that  the  maximum  number  of  modes  needed  to  represent 
the  far-field  at  the  angle  6^  i£  approximately  ka  sin  9. 

Applying  this  information  about  T^(y)  to  eq.   (40)  reveals  that  b^((j)^,z^)  can  be 
considered  the  double  Fourier  transform  of  a  band-limited  function  -  provided  the  separation 
distance  d  between  test  and  probe  antenna  is  large  enough  to  insure  the  values  of  R'^  do 
not  grow  extraordinarily  large  for  n^>(Ka)^  (see  eq.   (37)  or  eq.   (38)).     Consequently,  the 
sampling  theorem  [9]  can  be  utilized  to  convert  the  integration  on  the  right  side  of 
eq.   (40)  to  the  following  summation: 

^  .  -inct)-^  lYZ 

i; (Y)  =    I    b;(*J,  z^)e        °  e      °       ,  (42a) 

4kTk„aa      i , 
12  o 

1  ij 

where  (j)^  and  z^  refer  to  points  taken  at  increments  of  \^/2a  and  X2/2  radians,  respectively. 

4kTk„aa  i.l 
1  2  o 


Similarly,  I||(y)  be  written 


^  .  -in^-"  -lyz 

i;;(y)  =  — - —  I  K^'^i,         °  ^    °  •  (^2b) 


The  values  of  k^  =  2it/X^  and  k^  =  2t\I\^  need  not  be  chosen  much  larger  than  the  free- 
space  wave  number  k  =  Zit/X,  as  long  as  the  separation  distance  between  test  and  probe 

g 

antenna  remains  large  enough  to  keep  the        (y)  of  eq.    (37)  or  eq.   (38)  from  growing  extraordi- 
narily large  for  |n|   >   |Ka|.     It  is  difficult  to  get  a  reasonable  estimate  from  eq.  (37) 
or  eq.   (38)  of  how  large  the  separation  distance  d  should  be.     However,  experience  with 
planar  near-field  scanning  as  well  as  with  hypothetical  problems  which  are  solvable  analytically 
suggests  that  k^,  k2  can  be  chosen  approximately  equal  to  k  provided  the  probe  and  test 
antennas  remain  more  than  a  few  wavelengths  apart  (i.e.,  outside  the  reactive  fields)  as 
the  probe  scans  on  the  surface  of  the  cylinder  enclosing  the  test  antenna. 
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In  summary,  eqs.    (42)  show  that  the  amplitude  and  phase  of  the  output  of  the  probe (s) 
need  be  sampled  only  at  discrete  points  on  the  surface  of  the  scan  cylinder.     As  long  as 
the  probe  does  not  get  within  a  few  wavelengths  of  the  test  antenna,   the  data  point  spacing 
need  be  no  closer  than  about  X/2a  for  the  cj)^  angular  separation  and  X/2  for  the  incre- 
ments.    After  the  output  of  the  probe(s)  is  recorded  for  the  entire  cylindrical  grid  of 
data  points,  the  "fast  Fourier  transform"  (FFT)  can  be  utilized  to  compute  I^(y)  and  I^(y) 

from  eqs.   (42).     From  the  values  of  (I'(y),   I"(y))  and  the  receiving  functions  of  the 

s 

probe,  the  transmitting  function  T^(y)  of  the  test  antenna  can  be  evaluated  immediately 

from  eqs.   (27).     Finally  the  entire  electromagnetic  field  outside  the  test  antenna  is 

s  s 

determined  by  substituting  B^(y)  -  ^^^^^  ^o  eqs.    (4)  or  into  eqs.    (8)  for  just  the 

far-fields. 

6 .  SUMMARY 

The  basic  theory  and  formulas  needed  to  implement  the  probe-corrected  measurement  of 
antennas  by  scanning  on  a  cylinder  have  been  derived  by  a  systematic  approach  based  on  the 
cylindrical-wave  source  scattering-matrix  description  of  antennas.     This  approach,  which 
parallels  the  plane-wave  scattering  matrix  approach  used  by  Kerns   [3],  proves  simple  and 
straightforward  and  does  not  require  the  probe  or  test  antenna  to  be  reciprocal. 

The  derivation  began  by  expressing  the  EM  fields  radiating  from  the  test  antenna  in 

a  complete  set  of  cylindrical  vector  wave  solutions  (m,  N)  with  first  Hankel  functions 

(H^"*"^)  for  the  radial  dependence.     The  radiation  characteristics  of  the  test  antenna  are 

determined  by  the  weighting  coefficients  of  these  cylindrical  waves,  which  in  themselves 

remain  independent  of  the  particular  antenna.     Expressions  were  derived  for  the  far-fields 

(a  technique  was  introduced  in  appendix  A  for  deriving  the  far-fields  by  direct  manipulation 

of  the  Fourier  transforms  rather  than  by  the  more  involved  method  of  stationary  phase) , 

power  gain  function,  complex  polarization  ratio,  and  total  radiating  power  of  the  test 

g 

antenna  in  terms  of  the  cylindrical  waves  and  the  weighting  coefficients  B^(y) .  Thus 

the  problem  of  determining  the  radiation  characteristics  of  antennas  by  cylindrical  scanning 

reduces  to  that  of  determining  the  modal  coefficients  B^(y). 

If  an  ideal  dipole  antenna  were  available  to  measure  the  transverse  electric  field  on 

g 

a  cylinder  about  the  test  antenna,  the  coefficients  B^(y)  were  shown  to  emerge  immediately 
through  a  trivial  application  of  the  orthogonality  relationships  for  the  vector  wave  solutions 
(M,  N) .     When  an  arbitrary  probe  antenna  was  inserted  into  the  picture,  a  second  set  of 
cylindrical  waves  identical  to  the  first,  except  for  Bessel  functions  (J^)  replacing 
was  required  to  express  the  EM  fields  excited  by  the  sources  on  the  probe.     The  coefficients 

of  this  second  set  of  J    cylindrical  modes  along  with  the  input  to  the  test  antenna  were 

^  (1) 
related  to  the  coefficients  of  the  first  set  of  H        modes  and  the  output  of  the  test 

n  ^ 

antenna  through  the  "source  scattering  matrix"  for  the  test  antenna. 
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The  source  scattering  matrix  contains  a  complete  and  convenient  description  of  the 
transmitting,  receiving,  scattering,  and  reflection  properties  of  the  test  antenna.  The 

transmitting  elements  of  the  matrix  essentially  equals  the  required  modal  coefficients 

s 

B^(y).     Reciprocity  was  also  shown  to  display  itself  conveniently  as  a  relationship 
between  the  elements  of  the  scattering  matrix. 

A  similar  source  scattering  matrix  was  written  for  the  probe  antenna  but  with  reference 
to  the  origin  or  center  of  the  coordinate  system  situated  in  the  test  antenna.  This 
"common-center"  approach  along  with  the  use  of  ('^^^^^  j-^jj)  scattering  matrices  greatly 
facilitated  the  subsequent  derivation  of  the  transmission  formula,  which  contains  the 
required  coupling  product  between  the  transmitting  coefficients  of  the  test  antenna  and 
receiving  coefficients  of  the  probe  antenna.     Another  advantage  of  the  common-center  source 
scattering-matrix  approach  was  that  it  clearly  and  completely  separated  the  derivation  of 
the  transmission  formula  from  the  problem  of  determining  the  receiving  characteristics  of 
the  probe. 

The  characterization  of  the  probe  was  dealt  with  as  an  isolated  topic  in  section  4. 
It  was  found  that  the  receiving  coefficients  of  the  probe  could  be  either  measured  directly, 
or  transformed  from  its  measured  far-field.     The  specific  transformation  involved  depends 
on  the  choice  of  coordinate  system  (not  necessarily  cylindrical)  in  which  the  far-field 
of  the  probe  is  expressed.     To  conform  to  the  work  of  Leach  and  Paris   [2],  however, 
cylindrical  coordinates  were  chosen  to  express  the  far-field  of  the  probe  and  the 
transformation  was  accomplished  with  the  help  of  Graf's  addition  theorem. 

Finally,  it  was  shown  that  the  transmission  formula  could  be  inverted  to  yield  the 
transmitting  coefficients  of  the  test  antenna,  and  that  the  inversion  could  be  accomplished 
efficiently  through  the  use  of  the  sampling  theorem  and  FFT  algorithm.     For  reasonably 
low  Q  antennas  the  number  of  cylindrical  modes  needed  to  expand  the  nonreactive  field  was 
shown  to  depend  on  the  elevation  angle  (90-6),  and  was  given  approximately  as  ka  sin  9 
where  "a"  is  the  radius  of  the  smallest  cylinder  circumscribing  the  test  antenna. 
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APPENDIX  A.     DERIVATION  OF  THE  FAR-FIELDS  IN  TERMS 
OF  THE  CYLINDRICAL  WAVE  COEFFICIENTS 


In  this  appendix  the  far-fields  radiated  by  a  test  antenna  into  free  space  (eqs.  (8) 
of  the  main  text)  are  derived  in  terms  of  the  coefficients  of  the  cylindrical  wave  expansion 
(eqs,   (4)  of  the  main  text).     Although  Leach  and  Paris   [2]  have  derived  the  results  using 
the  method  of  steepest  descent,  the  alternative  derivation  shown  here  is  interesting  enough 
in  itself  to  warrant  a  separate  appendix. 

The  starting  point  for  the  derivation  is  eq.    (4a)  for  the  E  field  outside  the  test 
antenna  expressed  in  a  cylindrical  wave  expansion.     If  M  and  N  from  eqs.    (5)  are  substituted 
into  eq.   (4a),  and  the  large  argument  value  of  the  Hankel  function, 

.   T7 

,t(1)/     n  „    /    -nH     /2       i(Kp  -  -r) 

(kp)  %  (-i)  ^  h'  , 

n  /ttkp 

is  used  in  M  and  N,  the  far-electric-field  can  be  written 

iKp 


E(p,(J),z)  =  /  Ke  (p,(!))  e^'^^dy      ,  (Al) 

—00  ' 


where 


with 


IT  . 

.   —  — 2^  00 

S(P'^>  =  /i^^    '     I     (-)"[Bi(Y)M^B^(Y)M     ]e^"*     ,  (A2) 


M      =i(—  e  -Ke.) 
ny         p    p  if 

hi      =  h-rKe    -  ^  e,  +  K^e  )       .  (A3) 

ny  p      p      (j)  z' 

Symbolically,  eq.    (Al)  may  be  written 

E       =  F  ^[(KC  )  —  ]     ,  (A4) 
r^  y  K 

where  F  and          will  be  used  to  denote  the  Fourier  transform  and  its  inverse,  respectively. 
The  convolution  theorem  transforms  eq.    (A4)  to 

E(r)  =  F-^KC  ]F-lj^]d?  • 

j;->oo  — oo   ^  ' 

The  second  inverse  Fourier  transform  in  eq.  (A5)  is  nothing  more  than  an  integral  representation 
of  the  zeroth  order  Hankel  function  of  the  first  kind,  i.e. 

TT. 


_1    e'-^P        Te'-^P    iy(S-z)  (1)„  rr^ — — r^  /2i       4'"  ikr  -ik?  cos 

h-z^^^  =  ^  ^  ^  =  ^C^k/(5-z)^+p^)  ^X.y^  e        e  e 
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Thus,  the  far-electric-field  may  be  written 

I          —  — i  oo 

TT/— x        /I         4    ikr  t  r~l  r  ~  1  -ikC  cos  9  t  kc.\ 

^^^^  "  /2it^  ^       ^       Z  '^r  t^S^^  •  ^  ^ 

j;->oo  _oo   ^  ' 

But  this  last  integral  is  the  Fourier  transform  of  an  inverse  Fourier  transform.  Thus,  the 
result  is  the  original  function  (multiplied  by  211)  evaluated  at  k  cos  9,  i.e., 

/2Tf       4^  ikr. 


E(r)  =  e        e      k  sin  9  £  (Y=k  cos  6) 

(k  =  A^-Y""  =  k  sin  9  when  Y=k  cos  9)     .  (A7) 

Substitution  of  e^(Y=k  cos  e)from  its  definition  in  eq.    (A2)  gives  the  required  far-field 
expression  for  E  shown  in  eqs.   (8)  of  the  main  text, 

E(r,(}),9)  =  -2k  sin  9  ^ —    \  (-i)"[B^(k  cos  9)8  -iB^(k  cos  9)eg]e^"^  (A8a) 
e  x¥ 

H(r,c^,9)  =  .  (ABb) 
r^«'  o 


Finally  it  is  mentioned  that  the  same  technique  can  be  utilized  to  find  the  far-fields  of 
a  radiator  in  terms  of  its  plane-wave  spectrum,  i.e.,  to  derive  eqs.    (1.2-16)  of  reference 
[3b]. 
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APPENDIX  B.     DERIVATION  OF  THE  SOURCE  SCATTERING  MATRIX  EQUATIONS 

AND  RECIPROCITY  RELATIONS  FOR  THE  PROBE  AND  TEST  ANTENNAS 


Refer  first  to  the  test  antenna  of  fig.  1.     The  reference  surface  S     is  a  dividing 
plane  between  the  "test  antenna  proper"  and  the  shielded  "generator"  of  the  test  antenna. 
Let  Sp  be  the  surface  of  the  cylinder  of  radius  p  surrounding  the  test  antenna. 

In  the  region  between  the  generator  and  the  surface  S^  assume  that  the  EM  fields 
satisfy  the  following  very  general  form  of  Maxwell's  equations 

VxE  -  icoB"  =  0 
VxH  +  icoD  -  O'E  =  0 
D    =  e'E  +  T«H 

B    =  V'E  +  y"«H  .  (Bl) 

Consider  the  "adjoint  fields"  which  satisfy  the  adjoint  set  of  equations, 

VxE^  -  iwB^  =  0 

-:                                     VxH^  -  iojD^  -  a^-E^  =  0 

F    =  -  v^-H^ 

=  _7^.E^  +  7^.H^       ,  (B2) 

where  the  subscript  t  designates  the  transposed  dyadic.  A  straightforward  manipulation  of 
eqs.   (Bl)  and  (B2)  yields  an  "adjoint  reciprocity  lemma"  [3b] 


(i     (E^xH-ExH^)-nda  =  0     ,  (B3) 


S  +S 

o  p 

where  (E,H)  are  any  fields  satisfying  eqs.  (Bl)  and  (E^,  H^)  are  any  fields  satisfying 
eqs.  (B2). 

On  the  waveguide  reference  surface  S 

E    =  (a  +b  )e       ,     H^  =  n  (a  -b  )h^ 
_  ,  t  oo    o  t  oooo 

E^  =  (a^+b^)e      ,    If  =  n  (a^-b^)h 
t  oo    o  t  oooo 

f    e  xh  'nda  =  1  , 

S      °  ° 
o 


and  thus  eq.    (B3)  becomes 


2n  (b^a  -b  a^)  =  /     (E^xH-ExH^) 'e^da 
o     o  o     o  o  g  p 

p 
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(B4) 


Divide  the  fields  into  fields  emanating  from  sources  at  the  test  antenna  (denoted  by 
subscript  "o")  and  fields  emanating  from  the  probe  (denoted  by  subscript  "1").  Then 


/     (FxH-ExH^)-e  da  =  /     (FxH,-E  xH^+E^xH  -E,xH^)-e  da      .  (B5) 

g  P  gOlollolo  p 

p  p 

The  integrals  over  the  "o"  fields  alone  and  over  the  "1"  fields  alone  vanish  because  they 
satisfy  the  adjoint  reciprocity  lemma  separately.  Substitution  of  eq.  (B5)  into  eq.  (B4) 
gives  an  expression  which  will  be  called  the  source  adjoint  reciprocity  lemma, 

2n„(b^a  -b  a^)  =  /     (E^xH,-E  xH^+E^xH  -E,xH^) -e  da      .  (B6) 
ooooo        •'g      olollolo  p 

p 

If  the  E,  H  fields  in  the  coordinate  system  of  the  test  antenna  given  by  eqs.    (13)  are 
substituted  into  the  right  side  of  eq.    (B6),  it  becomes 

2        00  00  2 

%(^:VV:>=It-      I      /     (-l)M   [b^(Y)a-(-y)-b-(-Y)a^(Y)]dY  .  (B7) 
o        n=-<»    -oo  S=l 

To  get  eq.   (B7)  in  the  form  shown,  the  orthogonality  relations  (6)  between  the  M  and  N 
functions  and  the  Wronskian  for  the  Bessel  functions  were  used  ([7],  Formulas  9.1.15-9.1.17) 

Recall  that  ^^(Y)  3.nd  ^^^(Y)  coefficients  of  modes  with  arbitrary  sources  outside 

the  test  antenna  and  its  adjoint,  respectively.     Thus,  a^(Y)  and  a^^(Y)  may  be  chosen  inde- 
pendently and  arbitrarily.     Similarly,  a^  and  a^,  the  inputs  to  the  test  antenna  and  its 
adjoint,  may  be  chosen  independently  and  arbitrarily. 

In  particular,  if  we  choose 

a    =  1  a      =  0 

o  n 

eq.   (B7)  yields 

2      2        00  oo 

^■❖o+^^    I      I      I  (-l)V=(-Y)a^(Y)d.     .  (B8) 

O    O      S=l   n=-oo  -00 

But  this  is  nothing  more  than  the  source  scattering  matrix  equation  (14b)    for  the  test 
antenna  with 

r  =  b^ 

o  o 

„s  ,  .      4ttk:^     ,  , ,  n  ,  as  ,  , 
o  o 

In  a  similar  manner  the  second  source  scattering  matrix  equation  (14b)  emerges  if  we  choose 

a^  =  0    ,     a^^  =  6    6  6(y-6) 
o  n  sq  nm 
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In  other  words,  the  source  scattering  matrix  equations  (14)  derived  from  Maxwell' s  equations 
through  the  source  adjoint  reciprocity  lemma  (B6) . 

Equation  (B6)  also  reveals  the  reciprocity  relationships  between  the  source  scattering 
matrix  of  the  test  antenna  and  it  adjoint.     When  eqs.    (14)  are  inserted  into  eq.    (B7)  the 
necessary  and  sufficient  conditions  for  eq.    (B7)  to  be  satisfied  for  all  values  of  the  "a" 
coefficients  are 

=  r  (a) 

o  o 

R^^(_Y)  =  (_i)"  ^^t"(Y)  (b) 
-n"  nZkn' 
o  o 

T-<-,)  .  (-1)"  ^  (c, 


By  definition,  a  reciprocal  antenna  is  one  in  which  the  antenna  and  its  adjoint  are  identical. 
Thus  for  a  reciprocal  antenna  the  adjoint  superscript  "a"  in  eqs.    (B9)  are  removed  and  eq. 
(B9b)  becomes  identical  to  eq.   (B9c).     (Also,  eq.    (B9a)  becomes  a  trivial  identity.) 

We  have  just  completed  the  derivation  of  the  source  scattering  matrix  equations  (14) 
and  their  reciprocity  relationships  (15)  for  the  test  antenna  directly  from  Maxwell's 
equations  through  the  source  adjoint  reciprocity  lemma  (B6) .     For  the  probe  antenna  we  can 
proceed  in  the  same  fashion  as  for  the  test  antenna,  but  this  time  the  region  of  interest 

lies  between  S    and  the  probe    generator  which  is  separated  from  the  probe  antenna  proper 

p 

by  the  waveguide  reference  surface  S^.     Applying  the  adjoint  reciprocity  lemma  to  this 
external  region  yields  a  surface  integral  similar  to  eq.  (B3), 

I       (e'^xh'-  E'xH'^)-h'da'  =  0      ,  (BIO) 
J 

s'+s  +s 

o     p  oo 


where  the  primes  have  been  used  to  emphasize  the  fact  that  we  are  dealing  with  the  probe 
antenna.     Immediately,  the  radiation  condition  at  infinity  eradicates  the  integral  over 
S^,  and  eq.    (BIO)  becomes  identical  in  form  to  eq.    (B3) . 

Substitution  of  the  EM  fields  of  the  probe  waveguide  on        into  eq.    (BIO),  and  division 
of  the  fields  on        into  those  emanating  from  the  test  and  probe  antennas  yield  an  equation 
identical  in  form  (except  for  a  negative  sign)  to  eq.    (B6) ,  i.e.,  a  source  adjoint  reciprocity 
lemma  for  the  probe  antenna  is  derived.     Moreover,  substitution  of  the  E'  H'  fields  from  eq. 
(17)  gives  an  expression  identical  in  form  to  eq.    (B7)  but  with  the  primed  coefficients  of 
the  probe  system.     Continuing  the  derivation  as  for  the  test  antenna  gives  rise  to  the 
source  scattering  matrix  eqs.   (18)  for  the  probe,  and  reciprocity  relationships  for  the  probe 
identical  in  form  to  eqs.   (B9) .     That  is,  the  reciprocity  relationships  for  a  reciprocal 
probe  are  2^ 


^-n<-^>  =  ^'^'^^ 

o  o 

s:^^«^_^(-3.-Y)  =  s;^^(y,6)  .  (Bllb) 

Finally  it  is  mentioned  that  the  derivations  for  the  plane-wave  scattering  matrices  and 
their  reciprocity  relationships  for  antennas  have  been  extended  to  electroacoustic  trans- 
ducers [10].     The  analysis  in  reference  [10]  can  be  utilized  immediately  to  extend  the 
results  of  this  appendix  to  the  cylindrical  scanning  of  electroacoustic  transducers. 
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GLOSSARY  OF  SYMBOLS 


Radius  of  smallest  cylinder  circumscribing  the  test  antenna. 


a  ,a    or  a 
o    o  o 


b  ,b'  or  b" 

^  o    o  o 


n    n  In 

s     's  's 
b  ,b  ^,b,^J 
^  n    n  In-" 


Ingoing 
Outgoing 


*  wave-coefficients  in  the  test  and  probe  antenna  feeds,  respectively. 


Outside 
Inside 


source  cylindrical  mode  coefficients  for  the  test  antenna,  for 
'    the  probe  with  respect  to  the  C  coordinate  system  and  for  the 
probe  with  respect  to  the        coordinate  system,  respectively. 


Radiated  cylindrical  mode  coefficients  for  the  test  antenna. 


C(p  ' ,  (J) ' ,  z  ' ) ,  C^(p^,(})^,z^)   :     The  two  cylindrical  coordinate  systems  fixed  in  the  probe; 
the  first  (C)  centered  on  i 
on  the  probe  (see  fig.  3). 


the  first  (C)  centered  on  the  test  antenna  and  the  second  (C^)  centered 


d  : 


Distance  between  the  z-axes  of  the  C  and        system  (see  Fig.  3) 


(E,H) , (E' ,H' ) , (E^,H^) :     Complex  electric  and  magnetic  fields  referred  to  in  the  test  antenna, 
the  probe-C,  and  the  probe-C-j^  coordinate  systems,  respectively. 

e     :  Unit  vector  for  coordinate  a. 


Fourier  transform  parameter  for  the  z-part  of  the  cylindrical  modes 

(-00   <   Y   <   oo)  . 


r  ,r 

o  o 


n  n 


The  test  and  probe  antenna  input  reflection  coefficients  defined  at  the 

surface  S    and  S  ,  respectively, 
o  o 


Cylindrical  Hankel  functions  of  the  first  and  second  kind. 


,(1)' 


Derivative  with  respect  to  argument  of  the  Hankel  function  of  the  first  kind. 


I',   I"  : 

n  n 


Inversion  integrals  of  the  transmission  formula  (defined  by  the  right 
side  of  eq.    (26a)  and  (26b),  respectively). 


Cylindrical  Bessel  functions, 


Free-space  wave  number  =  2tt/X. 
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2    2  ]  /2 

<    :  -y  )        ,  taken  positive  real  when  y  <  k,  and  positive  imaginary  when 

Y  >  k. 

\     :  Free-space  wavelength. 

(M     ,  N    ),  ^it^^   '•     Cylindrical  vector  wave  solutions  with  J    and  H^"*"^  radial 

ny      ny         ny       ny  n  n 

dependence,  respectively. 

:  Characteristic  admittance  of  the  propagated  mode  in  the  test  antenna  feed. 

(4)^,2^)   :  Cylindrical  coordinates  describing  the  position  of  the  probe  as  it  scans 

the  test  antenna  (see  fig.  2). 

(p,i|),z),   (r,(}),e)   :     Cylindrical  and  spherical  coordinates,  respectively,  fixed  in  the 
test  antenna  (see  fig.  1). 

s       's  "s  's 

R  ,  R      or  R     ,  R,     :     Cylindrical  receiving  functions  for  the  test  antenna,  for  the  probe 
n      n  n  In 

with  respect  to  the  C  coordinate  system,  and  for  the  probe  with  respect  to 
the  C^  coordinate  system,  respectively. 

s     's  "s  's 

T  ,T      or  T     ,T,       :     Cylindrical  transmitting  functions  for  the  test  antenna,   for  the 
n    n  n  In 

probe  with  respect  to  the  C  coordinate  system,  and  for  the  probe  with 
respect  to  the  C-j^  coordinate  system,  respectively. 

gSq    g  sq     _       Scattering  functions  for  the  test  and  probe  antenna,  respectively, 
nm  nm 

I 

S^,  :       The  waveguide  reference  surfaces  for  the  test  and  probe  antenna,  respectively. 

OJ  :      Angular  velocity  in  the  suppressed  time  factor  e  ''"^'"(aj  >  0). 
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Figure  2.     Schematic  of  probe  antenna  and  its 
coordinate  system. 
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,rasJ-BYa  9:t,snxb'xooo 

Figure  3.     Schematic  of  the  two  cylindrical 
coordinate  systems  fixed  in  the 
probe. 
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.(O^'H  xsM  fi3c£) 
,1  ,a08  ,(.2.U)  .hnsJc 
\d'\L-^d\L  .qq  ,11  :J^b1 


to  liiiO'Ii}'. 


n  (order  of  Hankel  function) 


j±3syoDBo:t3.09is 


Figure  4.     Transmitting  coefficient  becomes 
negligible  not  far  outside  the 
circle  of  radius  ka. 
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